Mass Effect in Single Flavor Color Superconductivity 
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The nonzero strange quark mass effect in different types of single flavor color superconductivity 
and the phase diagram in a magnetic field are studied. We have obtained simple analytical forms of 
the quasi-particle energies for an arbitrary mass and explored the mass correction to the pressure and 
the transition temperature, ft is found that the mass reduces the pressure and transition temperature 
of strange quarks, but it doesn't change the ranking P n < Pa < fpoiar < fpianar < -Pcsl of the 
pressure for the four canonical single flavor phases. The phase diagram with magnetic field and 
temperature for a system of three flavors is obtained for two different values of the strange quark 
mass. The changes from the one obtained previously under the approximation of massless strange 
quarks are examined. 

PACS numbers: 12.38.Aw, 11.15.Ex, 24.85. +p 



INTRODUCTION 



Quark matter at sufficiently high baryon density and low temperature becomes a color superconductor(CSC) pj. 
CSC is characterized by a diquark condensate, which is analogous to the Cooper pair in an ordinary superconductor, 
but the structure of the condensate is much richer because quarks have the nonabelian color and flavor charges. 

The structure of the CSC states depends sensitively on the number of quark flavors and their masses [2-5]. For 
very high baryon density, where the masses of u, d and s quarks can be ignored, the ground state is in the color-flavor- 
locked(CFL) phase [6[, where quarks of different flavors pair. The situation becomes more complicated in moderate 
density because of the strange quark mass, /3 equilibrium and the charge neutrality conditions. A substantial Fermi 
momentum mismatch among different quark flavors is introduced and thereby reduces the available phase space for the 
cross-flavor pairing, such as CFL. Different exotic scenarios for cross-flavor pairing proposed in the literature (gapless 
CSC, LOFF CSC etc.) either run into various instabilities 043 or reduce significantly the condensation energy. This 
\ makes the single flavor pairing, which is free from the Fermi momentum mismatch and the instabilities, a competing 
ly*-) ■ alternative even though the pairing force here is expected to be weaker. There are a number of different paring states. 
The ones fre quently discussed in the literature include the spherical color-spin-lock(CSL) and nonspherical planar, 
polar and A [Ty-[l3 - Here, the adjective "spherical/nonspherical" refers to the symmetry of the order parameter 
under a space rotation. The CSL pairing is energetically most favored in the absence of a magnetic field, but the 
situation changed when a magnetic field is applied. The single flavor color superconductivity may be realized in the 
interior of a compact star during the later stage of its life, where a mag netic field is present. 

The presence of a magnetic field in the interior of a compact star [13| will offset the energy balance among the four 
canonical single flavor pairings. The spherical CSL phase has an electromagnetic Meissner effect [llj ]. but nonspherical 
phases: polar, A and planar phases do not. So if a quark matter of single flavor parings cools down through the critical 
temperature in a magnetic filed, forming CSL state will cost extra work to exclude magnetic fluxes from the bulk. 
Therefore, the magnetic contribution to the free energy may favor the nonspherical states. In a previous work, we 
have explored the consequences of the absence of the electromagnetic Meissner effect in a nonspherical CSC phase of 
single flavor pairing [lij and have obtained the phase diagram with respect to the magnetic field and the temperature. 
We found that under the plausible magnitude of the the magnetic field inside a compact star, the most favored state 
is not always CSL and nonspherical pairing states may show up. For the sake of simplicity, we considered both the 



infinitely massive limit and the massless limit strange quarks in 1AJ. The former limit is unrealistic given the typical 
chemical potential /i around 500MeV, the latter requires the mass of strange quarks, m s , to be much lower than quark 
chemical potential. On the other hand, m s has to be sufficiently large in order to win the competition with exotic 
cross- flavor pairings such as gapless CSC and LOFF. Both requirements may be compromised marginally for the value 
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of m s in vacuum (~ 150MeV) but will be problematic when the value of m s in medium becomes comparable with 
/i as was suggested by some numerical works such as [l5( . All these concerns warrant a systematic treatment of the 
single flavor pairing with an arbitrary quark mass. So we did in this paper. 

In the present work, we shall give a detailed investigation of the phase structure. For this purpose, we formulate 
the single flavor CSC for an nonzero quark mass in terms of the same NJL(Nambu-Jona-Lasinio)-like effective action 
employed in [3] and introduce the mean- field approximation for an arbitrary mass in section 2. Unlike the ultra- 
relativistic limit, where the cross-helicity(transverse) pairing dominates, the nonzero quark mass couples the cross 
helicity pairing channel and the equal-helicity(longitudinal) pairing channel and thereby complicates the gap matrix 
underlying the excitation spectrum. Fortunately, as will be shown in section 3, the gap matrix for an arbitrary mass 
can still be diagonalized analytically for all four canonical phases and our results interpolate both the ultra-relativistic 
limit and the non-relativistic limit in the literature. The ranking of the condensation energy in the massless limit 



remains intact when a nonzero quark mass is switched on. In section 4 we generalize our analysis in (14| of a three- 
flavor quark matter beyond the ultra-relativistic limit. Because the transition temperature of the nonzero m s strange 
quark paring is reduced, phase diagram with respect to temperature and magnetic field contains a region where only u 
and d flavors condensate. The size of this region is tiny for m s ~ f 50MeV but cannot be ignored for m s ~ fj,. Finally, 
we summarize our results and remark on some open issues in section 5. Throughout the paper, we shall assume zero 
masses for u and d quarks as we did in [13 . All gamma matrices are hermitian according to our notation. 

2. THE HAMILTONIAN UNDER MEAN FIELD APPROXIMATION 

In this section and the next one, we shall formulate the single flavor Copper pairing with an arbitrary quark mass. 
The Lagrangian density of the NJL-like effective action reads [l6| : 

C = $(-7„$, + m + - G^ v T l ^ v T 1 ^ (2.1) 

where T l = with X 1 the Zth Cell-Mann matrix, m is the quark mass and /i is the chemical potential. We set the 
effective coupling G > 0, in accordance with the interaction mediated by one-gluon exchange at high density and that 
mediated by instantons for intermediate density. The corresponding Hamiltonian is 

H = J d 3 r • V + m - fi^i/j + G^ v T l ^ v T 1 ^ (2.2) 

Like QCD Lagrangian, the diquark scattering in (|2.ip conserves the eigenvalues of 75 of each quark. At m = 0, the 
eigenvalue 75 coincides with the helicity so that the helicity of each quark is also conserved during the scattering. The 
process like 

{R,R)^{R,L) (2.3) 

with R(L) the right(left) hand helicity will never occur and the transverse pairing will not couple with the longitudinal 
one. For m 7^ 0, however, the helicity is not the eigenvalue of 75 and is no longer conserved. The two types of pairing 
do couple via (|2.3[) . 

The thermodynamic pressure 

P=^\nexp(-^) (2.4) 



n 1 \ t, 

with T the temperature and f2 the volume of the system and the ensemble average of the operator O is given by 

<0> ~ Tr[ex P (-f)] (2 ' 5) 

In terms of the plane- wave expansion: 

^ = ^ E( a P^P^ e * Pr + & P>P.* e ~ ipr ) ( 2 -6) 

with s(= ±5) the helicity defined by 

cr ■ pu PiS = 2spw p , s a ■ pwp, s = -2spw Pi5 . (2.7) 
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the interaction Hamiltonian reads 

1 

n 



= q 53 a ^s^' a -p.*2 a - p /, s ^^,si' lt p', s ' 1 7i/ u -p,s2 u lp/, s ^7i/ u p,si (2-8) 



The formulae 



p,p' 

the terms containing antiquark operators, 6's 



T ij T km = -^ijhm ~ S im S kj ) + -(SijSkrn + S lrn S kj ) (2.9) 

enables us to decompose the diquark interaction into color-antisymmetric and symmetric channels and the interaction 
within the former is attractive and therefore responsible for Cooper pairing for G > 0. We have 

a P '.s' 1 r ' a -p. ;i 2 a -p' !;j ^' a p,si = ^ a P ', s ' 1 eC fiL p ', s ^_ P!;52 £ c ap !Sl + the color symmetric interaction (2-10) 

where the 3x3 antisymmetric matrix e c in color space is defined by (e 1 ) = A5, (e 2 ) = A 7 and (e 3 ) = A2 and they 
coincide with the matrix representation of the angular momentum operators of spin one with respect to cartesian basis. 
We shall designate J x , J y and J z for A5, A7 and A2 below. Furthermore, the gap energy associated to spin-one Cooper 
pairing is expected much smaller than the chemical potential. Therefore, we may drop the antiquark contribution and 
keep only the color antisymmetric interaction under the mean field approximation. The relevant Hamiltonian takes 
the form 

ffeff = ^2v F (p- k F )a+ s a ps - — MA^uB2(p',P)%,^ i £ c a t _ pl y 2 a- P;S2 e c a PtSl (2.11) 

P> S P^P'^i,S2,Sl,S 2 

where 

^<,4; Sl , S2 (p',p) = M P ', s ' 1 747^'"-p. S 2'"- P ',4747^"p, ;il (2.12) 



and the approximation \J p 2 + m 2 — fi ~ vp(p — kp) has been made with the Fermi momentum kp = \J /i 2 — m 2 and 
the Fermi velocity vp — kp/fi. 

To simplify (|2.1ip further, we employ the explicit form of the four component spinor in the chiral representation 



E+2sp j 
IP ^rp,s 



2 p Vp,s 



where the two component spinor <j> given by 



cos § \ / -e"^ sin § 

e^sinfj C osf 



(2.13) 



(2.14) 



with (6, <p) the polar angles of p. Our choice of the phases of <fi p ± 1 is to make them corresponding to the two columns 
of the standard Wigner D-matrix of the angular momentum J = 1/2, i. e. 

(<t> p ,iA P ,-i)=Di(cp,d,-<p) (2.15) 

where 

Di >m {a, (3, 7) =< Jm'\e- u > a e- u ^e-' u >" i \ Jm > (2.16) 
with J the angular momentum operator and (a, j3, 7) Euler angles. The chiral representation of gamma matrices is 



<j v 
o v 

where a v — (1, a),a v = (1, — cr) with cr's the Pauli matrices. 



(2.17) 
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After some algebra detailed in the appendix A, we find that 



AG 



H cS = 5>F(P - k F )a+ s a ps - — £ '^(pOW 



p.s 



p,p' 



where 



$ » = E (-ir^e-^B slS2 (p) 



1 



1 1 

2 2 

-s 2 si s 2 - si 



(2.18) 



(2.19) 



with the phase # p . s is defined by the relation 



e- <9 "- = -*(-l) s -^ Pi _ s 0- P , s = ie 2 *^ (2.20) 

and Bi i(p) = B_ i _ i (p) = A = 22 ~ ^,Bi_i(p) = B_ n (p) = 1 . The repeated indexes in the second term of (|2.21j) 

are summed over with fj,, v — 0, ± and the summation ^ p extends to half of the momentum space. We have defined 

J± = (e 1 ± is 2 ) and J = £ 3 in (|2"38]l . 

Introducing a long range order (a_ pS2 e c a psl ) and expanding the interaction term of (|2.19[) to the linear order of 



the fluctuation a. 



pS2 e c a psl — (a_ pS2 e c a psl ), we obtain the linearized mean-field Hamiltonian 

h mf . = ^»f(? - ^Xa, + J^a;a; - 3^[a;$;( p ) + a;^( p )] 



p.s 



where the order parameter A v is defined by 



3A- 
4G 



(2.21) 



(2.22) 



and will be regarded the element of a 3 x 3 matrix with fJ,(v) the row(column) index. In terms of the Nambu-Gorkov 
basis 



A + 



—w 1 _ 
e p - 2 a_„ i e 

Pi 2 



1 ~ 

" 2 a 



P> 2 P, 



-4 r 



(2.23) 



the Hamiltonian (j2.21j) takes the form 



H MF . = — A**A^ + V(P - M + £ 'Kka, 



(2.24) 



where 



hp = 



v F (p — k F ) M 

-v F (p-k F ) 



and the 6x6 matrix M is defined by 

M = V3A^J„ x 
with x the direct product and A = — . We have 



v 2 F (p-k F ) 2 + MJVft 









-v F (p — k F ) 2 — M^M 



(2.25) 



(2.26) 



(2.27) 



where both MM^ and M^M have identical nonnegative eigenvalues, 8 2 with n = 1, 6 and the quasi-particle energy 
reads E^.n = \/v 2 F (j> — &f) 2 + Replacing the hamiltonian iJ of (|2 .4[) by the linearized one of ()2.21[) . we end up 
with the pressure under the mean-field approximation 



(2.28) 



p = - JLajt a; - 1 y: '(v F (p - k F ) - Ep . n ) + ^ ln f 1 + cx p 

In the massless limit, A = 0, the off-diagonal elements of the 2x2 matrix in (|2.26[) vanish and we are left with only 
the transverse pairing. 



5 



3. THE THERMODYNAMICS OF THE SPIN-1 COLOR SUPERCONDUCTIVITY 

The polar, A, planar and CSL are the four canonical phases mostly discussed in the literature of the spin-1 color 
superconductivity. Each of them corresponds to a particular diagonal form of the 3x3 matrix A c a introduced in the 
last section. The thermodynamics will be discussed in this section for an arbitrary quark mass. 

To gain more insight to the geometrical structure of these spin-1 phases, we introduce the following two sets of 
spherical basis 

e± = =F^(x±iy) e = z (3.1) 

and 

e± = T^Z(6±i<p) e = p (3.2) 

where 9, tp and p are the unit vectors in the directions of increasing 9, tp and p of the spherical coordinates of 
momentum p, given by 

p = (sin 9 cos tp, sin 9 sin tp, cos 9) 

9 = (cos 9 cos tp, cos 9 sin tp, — sin 9) (3-3) 
tp = {— sm tp, cos tp, 0) 

The extra phase factor e ±lv renders e± nonsingular at the north pole, 9 = 0. It is straightforward to verify that 

D 1 a0 (tp,0,-<p)=e* a -e (3.4) 

and the gap matrix takes the compact form 

M = V3Ai( tf a * (3.5) 

where the indexes a and f3 run over either the spherical basis (|3.ip or cartesian basis x, y and z. 

Now we are ready to introduce the four canonical spin-1 phases in terms of the circular basis (|3.ip . with respect to 
which 

11 , N 

J = pJ-e + H — — J + e_ + J eo (3.6) 

V 2 v 2 

with J± = J x ± iJy and Jo = J z - Each of the canonical phases corresponds to a particular form of the 3x3 matrix 
in (|3.5p fwith a labelling the rows and (3 the columns). We have 

A(polar) = Adiag.(0, 0, 1) (3.7a) 

1 \ / 

A(A) = A | or A(A) = A 1 | (3.7b) 
0/ \0 

A(planar) = Adiag.(l, 1, 0) (3.7c) 
v2 

A(CSL) = -pAdiag-fl, 1, 1) (3.7d) 



V3 



where A is the gap parameter to be determined. 
Correspondingly, the gap matrix 



M(polar) = J 3 - A ( ^ si f J f X ™ 6 . ) (3.8) 
Kl ' V 2 V JqAcos^ -Joe^smfl ' 
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M(A) = v^A 



2 9 



—Jo COS 

\J e iv sm8 -J e^sin 



|j e^ sin 6» 

2«V . 



(3.9) 



M (planar) = -A f | (J _ e ^ + J+ + e -i V) sin / - J_ e 2 ^ sin 2 § + J+ cos 2 § ) (3A0) 



Mcsl = ^a( '; 7 :/ ) CUD 



and 

The operators J7± and j7o inside Mosl are defined by 

J± =e;-J = ±e ± ^(J e ±iJ^) Jo = e* -J (3.12 



with Jq — 6 ■ 3 and J^, = <^ • J. They satisfy the same angular momentum algebra as J± and Jo in (13.6[) . 

Though the gap matrices (|3.8p . (I3.9p . (|3.10p and p. lip looks complicated, analytical expressions of the eigenvalues 
of MM* or MtM can be obtained for an arbitrary quark mass. Parametrizing the eigenvalues by A 2 / 2 (0), we find 
that 

(1/8) (VA 2 +8±A) 2 (d 4 =2), ^A 2 (d 4 = 2) for CSL phase 

f2( a\ _ J ( 3 /4) (cos 2 + 1 + A 2 sin 2 6») (rfj = 4), 0(e?i = 2) for planar phase , , 

; ^ -* ~~ » (3/2) (sin 2 6> + A 2 cos 2 6>) (dj = 4), 0(d 4 = 2) for polar phase { ' 



(3/4) (l ± v / A 2 sin 2 (9 + cos 2 6>) (d< = 2), 0(dj = 2) for A phase 



where the integer inside the parentheses following each expression indicates the degeneracy of each distinct eigenvalue. 
The details of the diagonalization is shown in Appendix B. The function f(9) is ^-dependent for the polar, A and 
planar phases and we shall refer to these phases as nonspherical. The CSL phase will be referred to as spherical 
because of the constancy of its f(9). 

Then the pressure corresponding to (|2.28[) becomes: 

p = ~W A " " S £ J {vf{p - kp) E ^ + ^ di ln f 1 + exp ("T^) ) (3 - 14) 



where E p> i — y/vp(p — kp) 2 + A 2 f 2 (9). Here the index i labels the distinct eigenvalues in each line of A3. 13)) with 
di the degeneracy. The summation over the entire momentum space is restored owing to the symmetry of /(0)'s 
under space inversion. Maximizing the pressure with respect to A, we obtain the gap equation T — 0, which 

determines the temperature dependence of the gap, A(T), up to the transition temperature. 
In terms of the parameter t — A ^ , the gap equation takes the form ln = j+x$ with 



M*)=Et / MsnOfm / dx 
4 1 Jo Jo 



^/x 2 +t 2 f 2 (9)[eV* 2 + t2 f?( e ) + 1] 



(3.15) 



It follows that 



t 

The condensation energy density of the CSC is given by 



T = =^e"5^ (3.16) 



P s -P n = Ps (t)^ (3.17) 



with s labeling different pairing states and Ao = Acsl(0) when m s = and we have 

^l 2 -±^ + h(t) + 2 g M 



Ps (t) = vpe-^^'i—^- + h(t) + 2^- a-p] (3.18) 



7 



with a — |(1) for nonspherical(spherical) phase and 

S(t) = Y)£/ MsmO / ^ln[ e -V^+* 2 /. 2 W + l] (3.19) 

The curves P(T) may be plotted parametrically according to (|3.16p and (|3.18|) without solving the gap equation for 
T > 0, as we did in [14]]. The transition temperature T A is determined by (|3.16[) in the limit t — > with A(0) the 
solution of the gap equation at T = 0. We find that 

T A = (^) (1 -i^^ } T c (3.20) 
Ao 

where X is a UV cutoff for \ p — kp | in the momentum integration and is assumed to satisfy the condition Ao << 
K « kp- We set K = 27MeV for the numerical calculation in this paper. For a given mass, T A is universal to 
all four phases and the ratio between the pressures of different phases is independent of the cutoff K . This cutoff 
matters only when we compare the gaps and pressures of different mass values. With Fermi velocity vf = vl — A 2 , 
T A is a monotonic decreasing function of A for < A < 1. Then the transition temperature with massive quarks 
is always lower than that in the massless limit. The factor p s vanishes at the transition temperature T C A . When 
A = 0, the corresponding curves of p s {T) is the same with what we got in [l4[ in the ultra- relativistic limit. We have 
Pcsl = f , Ppianar = 0.98, /0p O iar = 0.88 and pa = 0.65 at T — 0, in agreement with the values reported in [12]. In the 
non- relativistic limit, we get p po iar = Ppianar = 2<Pa = ^-5— Pcsl> consistent with the results in [Io| . 

The factor p s versus T/T A is plotted in Fig.l for A = 0.3 and A = 0.6. For p, = 500MeV, the former corresponds to 
to = 150MeV and the latter to m — 300MeV. The curves of Fig Q] implies the same inequality 

P n <P A < P P oi ar < P planar < Pcsl (3.21) 

as in the massless limit(A = 0). We expect (|3.21j) to hold within the whole domain of < A < 1. 







FIG. 1: The scaled condensation energy dependence on temperature with different masses A = 0.3 and A = 0.6. 
In what follows, we shall identify Acsl(0) with that of the one-gluon exchange [ToL Il7j. 

extrapolated to p = 500MeV and a s = J- = 1 with g the QCD running coupling constant. We can obtain the 
transition temperature T A=0 = ^-^Ao for u and d quarks in MeV. The transition temperature of s quarks, of T A ^° 
follows from (IX2TJ|) . For K = 27MeV, we hnd T c a3 = 0.98T C ° for m s = 150MeV and T c 6 = 0.683T C ° for m s = 300MeV. 
We should notice that the screening effect underlying the formulas (|3.22p comes from all three flavors in the massless 
limit. This inconsistency, however, will not affect our order of magnitude estimation. 
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4. THE PHASE DIAGRAM IN A MAGNETIC FIELD 



The physics of a quark matter in a magnetic field has received increasing attention because of the presence of a 
strong magnetic field in a compact star or during a noncentral collision of heavy-ions. The phase structure of 2SC in 
a magnetic field has been investigated in [IH, [l9| . Equation of state for the CFL phase in a magnetic field and its 
implications for compact star models have been studied in [20]. For an ultra-strong magnetic field B, the spacing of 
Landau levels becomes comparable or lar ger than the quark chemical potential, i.e. yeB > /i, the magnetic field will 
impact on the pairing dynamics of CFL [21| . For the typical value of u(=500MeV), this requires that B > 10 18 G, 
which may be implemented inside some magnetar. It was shown in [21] that an ultra-strong magnetic field may 
enhance the energy gap of the CFL for v 'eB >> /i and induce a magnetic moment of a Cooper pair. At a weaker 
magnetic field, \J~eB » \i, a de-Hass van-Alphen oscillation of the energy gap has been found [2311241. Alternatively 
a domain structure may be formed because of the chiral symmetry breaking and the axial anomaly [18j. For a spin-1 
CSC, in additional to above possibilities, a magnetic field may offset the balance between the CSL and nonspherical 
phases, producing a rich phase structure with respect to the temperature and the field shown in our previous work 
|14| . This mechanism will be further explored below taking into account the nonzero mass of strange quarks. 

The discussions of proceeding sections imply a nonzero order parameter 

$ =< ^ c r c A c v> > (4.1) 

in the coordinate space, where ip is the quark field, ipe — 72"0* is its charge conjugate, A c with c = 2, 5, 7 is an 
antisymmetric Gell-Mann matrices and T c is a 4 x 4 spinor matrix. We may choose r 5 = T 7 = for the polar and A 
phases, r 2 = for the planar phase but none of r c 's vanishes for CSL phase. Depending on the symmetry of (14.11) . 
their responses to an external magnetic field are quite different. 

For CSL phase, the diquark condensate (I4.1[) breaks the gauge symmetry SU(3) c xU(l) em completely. But the 
Meissner effect for a nonspherical condensate is incomplete, because it breaks the gauge symmetry partially. Among 
the residual gauge group, which leaves the diquark operator inside (|4.1[) unchanged, there exists a U(l) transformation, 
-0 — > e~i x&6 ~ lq ^^p with q the electric charge of ijj, 9 = —2^/3q4> for the polar and A phases and 9 — 4\/3q(j) for the 
planar phase. The corresponding gauge field, A^ is identified with the electromagnetic field in the condensate and is 
related to the electromagnetic field A and the 8-th component of the color field A s in the normal phase through a 
U(l) rotation 

A^ = A M cos 7 -A 8 sin 7 

V M = i4„ sin 7 + A* cos 7 (4.2) 

where the mixing angle 7 is given by tan7 po i ar A = 2^/3q(e/g) and tan7 P i anar = A\/Zq(e/g) for planar. The 2nd 
component of (|4.2I) V = because of the Meissner effect and thereby imposes a constraint inside a nonspherical CSC, 
A 8 = — A M tan7, which implies that: 

B 8 = -Btan 7 (4.3) 

with B = V x A. 

Expressing the gauge coupling ip^/^eqA^ + A 8 As/2)^ in terms of A^ and its orthogonal partner V M , we extract the 
electric charges with respect to A in color space, 

' '' = diag.(0,0, 1) for polar and A 



-p====diag.(l, 1, -1) for planar 
yjg 2 + 48<j 2 e 2 

Because of the nonzero charges of pairing quarks, the planar state is subject to the impact of Landau orbitals in a 
magnetic field, like that for CFL. 

The thermal equilibrium in a magnetic field Hz is determined by minimizing the Gibbs free energy density, 



-P+-B 2 



^{B'f-BH (4.5) 



2 

1=1 



with respect to A, B and B 1 . Ignoring the induced magnetization of quarks, the pressure P is given by (12. 4[) . with 
A given by the solution of the gap equation. Ignoring the induced magnetization of quarks, the pressure P is given 
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TABLE I: 





I 


ii 


III 


IV 


2 flavor 


CSL U , CSLd 


(polar) u , (planar) d 


(normal) u , (polar)d 


(normal) u , (normal) d 


3 flavor 


CSL U , CSLd.s 


(polar)u, (planar) d , s 


(normal) u , (polar) d , s 


(normal) u , (normal)d, s 



by (|2.4[) . with A given by the solution of the gap equation. For a nonspherical CSC pairing, the minimization with 
respect to B and B l is subject to the constraint (|4.3[) . For a hypothetical quark matter of one flavor only, we find 
that 

Gmm.j = -Pi - \vjH 2 (4.6) 

with j — n, CSL, polar, A and planar, where r\ n — 1, ?7csl = and ijj — cos 2 jj for a nonspherical CSC. The phase 
corresponding to minimum among £7 m i n 's above wins the competition and transition from one phase to another is first 
order below T c . 

The situation becomes more subtle when quarks of different flavors coexist even though pairing is within each flavor. 
Different electric charges of different quark flavors imply different mixing angles which may not be compactible with 
each other. Consider, for instance, a quark matter of u and d flavors with each flavor in a non-spherical CSC state 
with different mixing angles. (j4.3[) imposes two constraints, which are consistent with each other only if B = B s = 0. 
Then we end up with an effective Meissner shielding [IT| , making it fail to compete with the phase with both flavors in 
CSL states. On the other hand, one may relax the constraints by assuming that the basis underlying the condensate 
of u quarks differ from that underlying the condensate of d quarks by a color rotation. Consequently the constraint 
(I4.3[) for each flavor reads B s = —Stan 7" and B' s = —Bta,nj d . If both flavors stay in the polar or A phases, which 
allows B 1 3 to penetrate in, one may expect that an orthogonal transformation 

B' s = B s cos P - B 3 sin p 

B' 3 = B s sin /3 + B 3 cos p (4.7) 

could compromise both constraints. Such a transformation, however, lies outside the color SU (3) group and therefore, 
the mutual rotation of color basis is not an option. The phases of the two flavor quark matter (u, d) without Meissner 
effects, which can compete with (CSL, CSL), include (polar, planar), (polar (normal), normal(polar)), (A(normal), 
normal(A)) and (normal, normal). Notice the coincidence of the mixing angle of the polar state of u quarks and that 
of the planar state of d quarks. Also the normal phase does not impose any constraint on the gauge field and can 
coexist with any nonspherical CSC. 

The Gibbs free energies remain given by the equations of (|4.6p , but with P n and Pcsl referring to the total pressure 
of all quarks for normal and CSL phases. For nonspherical phases, P is the total pressure of all flavors with at 
least one of them in a nonspherical CSC state and 7 is their common mixing angle. For normal-CSC combination, 
7 refers to that of the CSC state. The number of combinations to be examined is reduced by two criteria: I) For 
two combinations of the same mixing angle, the one with higher pressure wins. 2) For two combinations of the same 
pressure, the one with smaller magnitude of the mixing angle wins. Because the function p s for various CSC phases 
also satisfy the inequalities (|3.2ip up to the transition temperature for an arbitrary mass, it follows that there are 
only four phases to be considered in each case of two and three flavors with nonzero quark masses, which are shown 
in Table I. 

The border between two phases are determined by the equation 

H 2 H 2 

P a + r, a — = P p + V0 — (4.8) 

with the subscripts a and p labelling the four phases I-IV. 

In a multiflavor quark matter the Fermi momentum of each flavor is displayed from each other to meet the charge 
neutrality requirement (The color neutrality condition is ignored owing to the small energy gap associated to the 
single flavor pairing). In what follows, we shall consider the quark matter of two massless flavors (m„ = rrid = 0) 
and a massive flavor (m s ^ 0), coexists with electrons. Within the mean- field approximation employed in proceeding 
sections, the Fermi-momentum displacement can be determined in the ideal gas limit at zero temperature. The total 
pressure under this approximation reads 

/ / v / / 
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where Ef, n/ and are the kinetic energy density and number density of the quark flavor / with / = u,d,s and 
qf — (2/3,-1/3,-1/3), E e and n e are corresponding quantities for electrons. A charge chemical potential (i q is 
introduced with the (...) of (|4.9p the charge number density. We have 



q rkf l4 

4^ 



£ / = ^/ dppzjpi+mj E e = -^ (4.10) 



n/ = ifc| ne = i fc ' (4 -"» 



with mf = (0,0, m s ). The Fermi momenta, fc/ and k e are determined by the equilibrium conditions 



\ dk f J ~ V dk e J , , 

and the charge neutrality constraint 



= (4.12) 



^g /n/ -n e = (4.13) 
/ 

We find that k u = 1.001/x, fc d = 1.01 and k s = 0.941^i for m s = 0.3/x, and that k u = 1.004/x, fc rf = 1.039/i and 
fe s = 0.744/i for m s = 0.6/i. We got this H-T diagram, Fig|3J and H is defined by 

Ho = ^ (4.14) 




FIG. 2: H-T phase diagram. These four diagrams corresponding to the ultra-relativistic limit m s = 0, m s » fi and m 3 = 
150MeV, m s = 300MeF. 
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When the strange quark mass is ignored, the transition temperature of all flavors are the same, this three flavor 
phase is what we got in the diagram "A = 0". The transition from one to another is first order. When we assume 
a large mass of s quarks, m a >> /i, This two flavor case is certainly unrealistic. When m a is comparable to //, the 
transition temperature of strange quark parings is reduced, so when T C A < T < T£, strange quarks become unpaired, 
while u and d quark parings remain. Below T C A , the condensation energy of strange quarks rises like (T C A — T) 2 . The 
transition from three flavor CSC to two flavor CSC is therefore of second order at zero magnetic field. Since only the 
condensation energy, not its derivatives, enters (14.81) . both the phase boundaries and their slopes with respect to the 
temperature are continuous at T C A . The dashed line in Figgis the border between the three flavor region and the 
two flavor region in the phase diagram. In the left region of the dashed line, the corresponding phases combination 
is "3 flavor" in the TableJH and the right region is corresponding to the upper line "2 flavor". The values of the 
bounder between two regions is close to these in massless limit in three flavor, also these nonspherical phases occupy 
a significant portion of the H-T phase diagram for a magnitude of the magnetic field of order 10 14 G. This strong 
magnetic field is plausible in a compact star. Due to the not apparent change to the border between the possible 
phases, the physical implications of the mass effect to the cooling behaviors or the latent heat released as the star 
cools through the phase boundaries won't be discussed in this paper, it won't give corrections in order with [22]. The 
smallness of the spin-1 gap makes CSC being of type I and the critical magnetic field B ~ /xAo << /i 2 /e. Therefore, 
the magnetic impact on the pai ring dy namics as well as the quark matter magnetization may be neglected [22j , unlike 
the situation considered in [2l"l. [23 . [24 1 . 



5. CONCLUDING REMARKS 



For the quark matter of moderate density, such as that may reside in the core of a compact star, the mass of strange 
quarks, m s , is not much smaller than the quark chemical potential fi and need to be taken into account. In this paper, 
we extend the study of the four single flavor phases of color superconductivity to include the effect of the nonzero 
m s . In spite of the complication of the coupling between the cross-helicity and equal helicity channels, the excitation 
spectrum is obtained analytically. We have explored its correction to the pressure and the transition temperature 
numerically. It is found that mass effect reduces the pressure and transition temperature of strange quarks, but it 
doesn't change the ranking P n < Pa < -Ppoiar < fpianar < -Pcsl of the pressure for the four canonical single flavor 
phases for the values of m s we examined. We suspect that the above inequality holds for an arbitrary m s . Then we 
generalized the previous work to the quark matter with massless u and d flavors and the massive s flavor. Because the 
transition temperature for strange quark parings is lower than that of massless quark pairings, this new H-T diagram 
consists of the two flavor CSC for T C A < T < T c ° and the three flavor CSC for < T < T C A . The three flavor and 
the two flavor region will be occupied by the same phases I- IV in Table I with the same relative positions. There is 
a second order phase transition from three flavor condensate to two flavor condensate at transition temperature T C A 
where strange quarks condense. The phase boundaries in the two regions join smoothly. As an order of magnitude 
estimate, we calibrated our model T c against the that from the QCD one-gluon exchange in the chiral limit and found 
the typical magnitude of the magnetic field in the phase diagram falls within the range of the plausible magnetic field 
inside a compact star in the literature. 

On the other hand, the effective Lagrangian (|2.1[) we employed in this paper is by no means the most general one. 
The Lorentz covariance of (|2.1[) is unlikely in a medium and the coupling G may depends on m s . Taking the one-gluon 
exchange as a reference, the screenings of the color magnetic channel and the color electric one by the medium are 
very different and should depends the masses of quarks. These properties is likely to persists qualitatively at the 
moderate density and should be reflected in the effective action to some extents. Also, the inequality (|3 . 2 1 1) may not 
be as robust as people thought. A purely Ginzburg-Landau analysis (2f| reveals some parameter region where the 
ranking (13.211) is offset even without a magnetic field. The microscopic mechanism supporting this observation, based 
on the most general four-fermion effective action or others, remains to be unveiled. 

Throughout this paper, we take the massless limit of u and d flavors and then gapless excitations exist in all I- IV 
phases of Table. I. In reality, the chiral restoration transition from low density to high density may be a crossover and u 
and d quarks may also acquire nonzero masses from the chiral condensate. Consequently, the excitations of the phase 
I where all flavors are in CSL will be gapped. This is welcome since it will slow down the direct Urea processes of 
cooling in a compact star by spin-1 CSC alone [16j |. But the gapless modes remains for nonspherical states and phase 
diagrams Fig. 2 are still valid qualitatively. Therefore the magnetic field inside the star cannot exceed the magnitude 
along the border line between I and II of Fig. 2 for a slow cooling process. 
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Appendix A 

In this appendix, we shall fill in the details from (|2~TT|) to ([2718]) . Substituting ([2713]) and (f2TTT|) into (|2~T2|) . we 

obtain that 

A s ' ltS ' 2 ; S1 , S 2 (P 1 , P) = 2CV iS ^ (p')C_ S2 _ Sl (p)^, ;S , a u (f>- PtS2 (l) i _ p/ ^a v (t) PtSl 

+ 2C_ si _ 4 (p')C S2S1 (p)^ /s , ( 7^_ p , S2 t _ p , iS ,a^ PiS1 (A.l) 
+ 2C 4 _ 4 (p')C_ S2Sl (p)0 p , !S ,a>_ p , S2 t _ p , s ,a^ p , Sl 
+ 2C_ sl4 (p')C S2 _ sl (p)0 p , iS ,<T ,/ 0_ P)S2 t _ p , iS ,r</»p, Sl 



where C ss i (p) = ^ ' ^ E+2s ^ E+2s p ) follows from the identity 



(C i7 ')a^(cTj) 7 5 = 2<5 Q< 5<5 ( 3 7 — SapS-yS (A-2) 



that 



b l',s' 1 dV< t , -v,s 2 4>-p.^ v ^P,s 1 = tfpi^o 1 



'<f>P,s a ^- p > , s -y '4>- 


P,Sl 


, -P,» 3 ^ip/y a ^P,« 1 




^-P^-p',^^ 


'P: S 1 


'p.si^Lp'.s^-p,^ 





(A.3) 



For two different momenta, p and p', we have 

<f>l,s'<t>r,s = (D^(^,6i,- v i)dHv,9,-^) s , s = D], s (R) (A.4) 

where R stands for the Euler angles corresponding to the product of the rotations specified by (tp, 9, —ip) and 
(<// , 9', —ip'). Together with the orthonormal relation cf> P S P , S ' = S ss r, we obtain that 

A s > 1 , s > 2 :s 1 ,s 2 (p'j P) 

+ 4(-l)^+^- 1 C s; _ 4 (pOC Sl _ S2 (p)e l( ^ p - + ^S^J,_ S2 ( J R)^ 4si (i ? ) (A.5) 

where the phase 6* PjS is defined in (|2.20p and satisfies the relation 

e i0- P ,s = _ e ie P<3 (A. 6) 

Because 

Dl_JR)Dl s ,JR) - D], si (R)Dl,_ S2 (R)) = detDi(R)e s ,_ s ,e sl - S2 = e s ,_ s ,e 8l . S2 (A.7) 

in (|A.5p . e s ' i _ s ^e Sl _ S2 ^ requires that s[ = s' 2 ,si — S2- Then the diquark operator of equal helicity is even in p 
because of the equation (|2.20[) . so sum over p will make it vanish on account of (|A.6[) . So this part doesn't contribute. 
Using the formula of Wigner D-functions 

A B G \ ( A B C\ 
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we find: 

11 i \ / l l 



D !> s ( R ) D \>s ( R ) = 3 ( 2 < 2 ' > ± ' ) ( 2 2 1 ) i R ) 

s i~ s 2 y ' -« 2 s i v ' V s\ —s 2 s^ — s^J V — s 2 si s 2 — Si y s 2 ''l' 32 51 v ' 



2 2 u 12 2 

Sj — S 2 S 2 — s'l / \ — S2 Si S 2 - S 



" )D^ S2 _ Si (R) (A.9) 



where Si , sL s\, s 2 can take values of ±i. The D°f , _ _ , (i?) part doesn't contribute either, because it also requires 
= s' 2 , s l = s 2j which makes it vanish when to sum over p. 
It follows from the second equality of (|A.4|) that 

Dl im2 (R) =Y, D ™rn 1 (<P'> e '>-<P') Dl m rn 2 (<P> e >-<P) (A.10) 
m 

and we arrive at 

^i, s - sl , S2 (p',p) - 12(-l)^+^- 1 C s; _ s ,(p')a i - S2 (p)e I 

X J2 D ks' 2 -s' 1 ( R ) Dl m,s 2 -s 1 (R) + - (A.11) 



! -S' 2 S 2 - S 1 J \ -82 Si S 2 - Si 



with the terms not contributing to the summation over momenta. 
Substituting (| A. 1 1|) into ([2. lip , the interaction term of (|2.11j) becomes 

E A 4i4 . SliS2 (p',p)4 /s , £ c aL p , !4 a_ p>S2£ c a PjSl = 12^S^(p')H;(p) (A.12) 
p,p',Si,s 2 ,si,s 2 p,p' 

with 

S » = E (-l) S2 ^e-^C Sl _ S2 (p) ( I ^ * ^ ) D^: 2 _ Si (^,e,-^)a_ pS2 J^ psi (A.13) 
If follows from the explicit form of the phase factor (|2.20[) and the symmetry properties of the D-functions that 

S^(-p)=E(- 1 ) S2 "' e " i9pS2C '- 1 ^(P)f 4 I )nl; s2 _ sl (<p,9,-cp)a- pS2 J"a psi (A.14) 



On writing J^ p , S^(p) with the summation ^ extending half space of and we end up with (|2.18[) with 

$-(p) = S»+H^(-p) (A.15) 

given by (l2~Tg|) . 



Appendix B 

In this section, we will give the details of the diagonalization procedure of the 6x6 matrix MM** for each single 
flavor phase. We shall write MM^ = A 2 A4. The eigenvalues of M. corresponds to f 2 (6) shown in (I3.13[) . 
The polar phase: 

It's straightforward to show 

_ 3 /(sin 2 # + A 2 cos 2 #)J 2 \ ( 

Mp^-^y (sin 2 9 + X 2 cos 2 9) J 2 ) [ti - L) 

and the color operator Jq decouples. The eigenvalues of Jq are 1, 1, and the functional forms of f(0) are therefore 
given by the 3rd line of (|3.13[) . 
The A phase: 



14 



In this case we have 

_ / (2sin 4 f + iA 2 sin 2 0)J o 2 -AsinSW 2 \ 

y -Asin6W 2 (2cos 4 f + iA 2 sin 2 0)J 2 J ^'^ 

The color operator J5 2 , which has eigenvalues 1,1 and 0, decouples again. The forms of f(9) given by the fourth line 
of (|3.13p correspond to the eigenvalues of the 2x2 matrix obtained from (|B.2j) by replacing Jq with its eigenvalues. 
The planar phase: 

The diagonalization of MM* is less straightforward because of the coupling between the helicity and the color 
indexes. In terms of J± = J±e f%,p , we have 

AW = l( a c b d ) (B.3) 

where: 



cos 4 B - + iA 2 sin 2 $\ J'_J+ + (sin 4 °- + -^A 2 sin 2 oj J+f_ - ±-(1 - \ 2 ){J n + J+) sin 2 (B.4a) 
b = -sm8e- iv [f + , J'_] c = 6 t (B.4b) 



2 



d = (^in 4 1 + I A 2 sin 2 |j JLJ'+ + (cos 4 6 - + h 2 sin 2 0) f+JL - \{l - A 2 )(Ji 2 + J?) sin 2 (B.4c) 

Since Jj_ = J-|_eT 1 ^ and Jo satisfy the same angular momentum algebra as J± and Jo do, we shall work in the 
representation where 

_ / 

J'_ = \/2 1 ) . (B.5) 
\° 1 

and Jo = diag(l, 0,-1). It follows that 

2sin 4 f + ^ sin 2 6> ±(1-A 2 )sin 

1 + cos 2 9 + A 2 sin 2 9 ) (B.6a) 

i(l-A 2 )sin 2 6» 2cos 4 f + ^sin 2 

1 

b = Xsm9e~ tv | ) (B.6b) 
0-1 




2 cos 4 | + %- sin 2 I 






1 + cos 2 9 + A 2 sin 2 i 



i(l-A 2 )sin 2 







2 sin 4 I + sin 2 I 



|(1 - A 2 )sin 2 6» 
By permutations of the rows and columns, this 6 by 6 matrix is transformed into the block-diagonal form 



(B.6c) 



M 



Planar 



1 + cos 2 9 + A 2 sin 2 I 







1 + cos 2 9 + A 2 sin 2 9 
M 4 



(B.7) 



where M4 is a 4 by 4 matrix, given by 



/ 2 sin 4 § + £■ sin 2 1 



Ma = 



±(1 - A 2 )sin 2 6> XsmOe-^ 



i(l- A 2 )sin 2 2 cos 4 
Xsm9e zv 




f + lfSin 2 f 



2 cos 





4 9 1 A 2 





-Asinfle-^ 
sin 2 6» |(1 - A 2 )sin 2 ( 



|(1 - A 2 ) sin 6* 2 sin 4 f + sin 2 0/ 



(B. 



15 



It is straightforward to show that the secular equation 

det(M 4 - z) = z 2 (z - 1 - cos 2 9 - A 2 sin 2 9) 2 

which, together with (|B.7[) . yield the eigenvalues in the second line of f|3. 13[) . 

The CSL phase: 

In terms of the operator J± and j7o, the matrix M. of CSL takes the form 

u _ 1 (J-J+ + \ 2 Ji \e-*[J ,J-]' 

M CSL - - ^ J +] J + J_ + X 2j2 

The operators J± and J§ satisfy the same algebraic relations as J± and Jq, such as [j7o,J7±] 
representation where J7q is diagonal, i.e. j7o = diag.(l, 0, —1), 



±J±. 



10 
J+ = y/2 [ 1 




and 




we have = diag(0, 2, 2), J + J_ = diag(2, 2, 0). Substituting these into (|RT0t . we find that 



■Mcsl = 



/A 2 





V o 




2 


V2Ae 








2 + A 2 








V2Ae-' 


2 + A 2 





\ 





-y/2e-^ 



2 

A 2 / 



By permutations of the rows and columns, this 6 by 6 matrix is transformed into the block-diagonal form 



/A 2 
OA 2 
2 

-V2Ae^ 

\ 






2 + A 2 










2 + A 2 -V2Xe-^ 
-V2Ae^ 2 



(B.9) 



(B.10) 
In the 

(B.ll) 



(B.12) 



(B.13) 



and the eigenvalues in the hrst line of A3. 13)) follow then. 
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